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Abstract 



The Balitsky-Fadin-Kuraev-Lipatov equation in the next-to-leading logarithmic ap- 
proximation is solved using an iterative method. We derive the solution for forward 
scattering with all conformal spins. A discussion of the infrared finiteness of the 
results is included. 



1 Introduction 



The Balitsky-Fadin-Kuraev-Lipatov (BFKL) pQ formalism enables the resummation 
of logarithms appearing in scattering amplitudes, which are large in the Regge limit, 
where the center of mass energy yfs is large and the momentum transfer y—t fixed. 
In this approach the high energy cross-section for the process A + B — > A' + B' is 
factorised as 

= / § / S W W ' ( k - k " A = ln i) ■ « 

where $a,b are process-dependent impact factors and / (k a , k&, A) is the process- 
independent gluon Green's function describing the interaction between two Reggeised 
gluons exchanged in the t-channel with transverse momenta k a ,b- m this Letter we 
use the Regge scale sq = |k a | |k&|, a different choice would modify the impact factors 
in such a way that the prediction for the cross-section remains unchanged. 

The representation in is valid [2 in the leading logarithmic approximation 
(LLA), where terms of the form (a s A) n are resummed, and in the next-to-leading 
logarithmic approximation (NLLA) [3], where contributions of the type a s (a s A) n are 
also taken into account. This formalism is valid in both the forward and non-forward 
cases j3]. In this Letter we deal with the former but the proposed method is also 
applicable to the latter. 

In recent years there have been many studies of the behaviour of the gluon Green's 
function in the NLLA 0. This Green's function is obtained as the solution of an 
integral equation where radiative corrections enter through its kernel. The main 
difficulty in solving the equation analytically in the NLLA stems from the logarithmic 
dependence introduced by the running of the coupling. In this work we show how it 
is possible to solve the equation using an iterative method directly in energy space, 
considering the full kernel with scale invariant and running coupling terms. In this 
approach we keep all the angular information from the BFKL evolution, solving the 
equation for a general conformal spin without relying on any asymptotic expansion. 

In Section El we present the BFKL equation in the NLLA in dimensional regu- 
larisation and show how, introducing a cut-off in the real emissions, it is possible to 
write it in a form suitable for iteration. In Section |3] we solve the equation using an 
iterative method and compare with previous results in the literature. In Section 0] we 
present our conclusions. 
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2 The BFKL equation in the NLLA 



To write the BFKL equation for the gluon Green's function, / (k , k&, A), it is con- 
venient to introduce the Mellin transform, f u (k a ,kb), in A space as follows 

1 ra+ico 

/(k ,k 6) A) = — / due^f^KM)- (2) 

ZIYI Ja—ioo 

In this way the BFKL equation in the NLLA for forward scattering can be written 
in dimensional regularisation (D = 4 + 2e) as 

uf u (k , k 6 ) = 6^ (k a - k 6 ) + J d 2+2 *k' K (k a , k') f u (k', k 6 ) . (3) 

The kernel in the NLLA is expressed in terms of the gluon Regge trajectory 
u>( e ' (k 2 ), and the real emission component 7j, K, r (k a , k), in such a way that 

K (k a , k) = 2 w « (k 2 ) 5^ (k a - k) + JCr (ka, k) . (4) 

It is convenient to split the kernel JC r into two parts: a e-dependent, K,\ e \ and a 
e-independent, K, r . In the integral over real emission we also perform a shift in the 
variable of integration over transverse degrees of freedom of the form k = k' — k a to 
write 

loU (k a ,k 6 ) = 5^ (k a -k 6 ) + J d 2+2 *k2u^ (k 2 ) 5^ (k a - k) U (k,k 6 ) (5) 
+ J d 2+2 *k (k a , k a + k) f u (k a + k, k b ) + J d 2+2 *k t r (k a , k a + k) f u (k a + k, kfc) . 

In order to explicitly show the cancellation of the poles at e = we split the 
integral over transverse phase space for K,fi into two regions separated by a cut-off 
A, i.e. 

uf u (k a , k b ) = 5^ (k a - k b ) + J d 2+2 ^k 2 (k 2 ) 5^ (k a - k) f u (k, k 6 ) 

+ J rf 2+2e k/c^ (k a , k a + k) (o (k 2 - a 2 ) + e (a 2 - k 2 )) f u (k a + k, 

+ J d 2+2 *kK r (k a , k a + k) f u (k a + k, k b ) . (6) 

For any finite A we now introduce an arbitrary dependence of O {jj^j by using the 
approximation f u (k a + k, kb) ~ f u (k a , k b ) for |k| < A. This dependence is negligible 
when the external scale k 2 is large. It is possible to introduce the cut-off in several 
different ways but all the choices must have the same A — > limit. In this way we 
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can approximate Eq. (JBJ) by 

c^(k a ,k 6 ) = 5^(k a -k b ) (7) 
+ {2 w « (k 2 a ) + J d 2+2 *k (k a , k a + k) 9 (A 2 - k 2 ) } f u (k a , k 6 ) 

+ J rf 2+2e k [Kf (k a , k a + k) (k 2 - A 2 ) + t r (k a , k a + k)} f u (k a + k, k 6 ) . 
In dimensional regularisation the gluon Regge trajectory can be written as [2j 



2CU& ( q 2 ) 



a 2 r 2 (l-e) f/? /l 



+ 



(4tt) 
'4 7T 2 5 /3b 



8 (40 2e [TV, 
• 2 ^ 32 



28 /3„ 



T + iFM7 + 2 "^)-T + 2C(3)-^, 



where /3q = y^c — 
in the MS scheme. 



Q'.s 
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9 iV, 



„2 

and a s (/i) = //is the renormalisation scale 



The e-dependent part of the real emission kernel reads [3] 



^ (q,q + k) 



a s fi 



-2e 



7T 



1+e 



1 f | a, r(l-e) 



(47r) e k 2 \ 4 (4tt 



' A 7T 



5/3 



32 



' /k 2N 
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28 ft 



1 - e z 
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9 N r 



(9) 



For our purposes we are interested in the integration of this piece of the kernel over 
the phase space limited by the cut-off, i.e. 



J d 2+2e kIC^ (q,q + k)#(A 2 -k 2 



(10) 



1 <*s_i /a 2 v r | a a r(i- 



r(l + e) (47r) e e \fi 



1 / A 2 N 



' I TV 2 5Po 

3 3 ' 3iV c 
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2\ e 
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When this term is combined with the trajectory, the poles in e cancel and we obtain 
a finite expression depending on A: 



co 



(q 2 ,A 2 ) = lhn|2^(q 2 )+|rf 2+2e k^(q,q + k)^(A 2 -k 2 )} 

1 rr 2 .>', 



(11) 



-« s ln- + T 



■ + 



3 



ln^-6C(3) 



To simplify our formulae we introduce the notation 

t 

A 2 



Wo (q 2 ,A 2 ) = -e(|q|A)ln^ + r ? 
3 



(12) 



with 



£(X) = a s + ^ 



5ft _ ft X 
3 3 3JV C N c n /i 2 



and 



(13) 



^ = ^oC(3). 



(14) 



In this way we can write 

lim J d 2+2e k (k a , k a + k) 9 (k 2 - A 2 ) f u (k a + k, k b ) = 

I d " k ^ ( k ') ° ( k ' " A ') fw (ka + k ' k& 



(15) 



With these conventions, the forward BFKL equation in the NLLA can finally be 
written as 



(o;-^ (k 2 ,A 2 ))^(k a ,k fe ) = 5^(k a -k b ) 



(16) 



d 2 k ( -j-^ (k 2 ) 9 (k 2 - A 2 ) + Kr (k a , k a + k) ) /„ (ka + k, k b 



where jHj 
IC r (q, q') 



3+1 



47T 



iq- q 
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^2 
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(17) 



X (q 2 + q' 2 + (q' 2 q 2 ) ^ q' 2 / 
j4\ \ \ 
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2^/2 
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q 

q'- 
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V 1 



1 — X 
ln(l - t) 



dt 



t 



- 1 



(q 2 - q' 2 ) 2 



(q-q') 2 (q + q') 2 , 



1 roc 



dz- 



lll 



q' — zq) 2 q 



/2 



It is important to note that this kernel contains the full angular information in the 
BFKL evolution. Writing the equation as in ()16j) is very natural in the sense that 
there is a clear separation between the virtual contributions on the left hand side, 
and the real emissions, integrated over phase space, on the right hand side. 



4 



To study the dependence on A of Eq.(16) we take the derivative with respect to 
A 2 of the A-dependent terms, i.e. 

^ {^o (K, A 2 ) f u (k a , k 6 ) + J d 2 k-^£ (k 2 ) 6 (k 2 - A 2 ) f u (k a + k, k 6 )} = 

^ (A 2 ) f u (k a , k 6 ) - -L d9t (A 2 ) f u (k a + A, k 6 ) . (18) 

For a sufficiently smooth f w (k a ,k&) this expression is small in the A — > limit. In 
fact, the approximation made in Eq.(7) is only good for a smooth f w (k a , k b ), and the 
A-dependence can ultimately be studied numerically. 



3 Iterative solution in the NLLA 

The BFKL equation in the NLLA for the forward written in Eq. (|16J) can be 

solved using an iterative procedure similar to the one applied in [H] for the LLA. We 
will, in the following, show how this works in the NLLA. In the NLLA it becomes 
meaningful to study the dependence on the renormalisation scale. Since there are in 
principle many different physical scales in the problem, there are many possible choices 
also for the renormalisation scale. Here we will assume that the renormalisation scale 
is chosen to depend on the arguments of the Green's function only 

/x = ^(k 2 ,k 2 ). (19) 

Other choices are possible and can be studied in this formalism. The /i-dependence 
of the trajectory and the kernel will be explicitly shown in all of our expressions from 
now on. It is also convenient to use the following notation for the kernel 

£ r (k a ,k a + k,A 2 , / u(k 2 ,k 2 )) = -^e(k 2 )^(k 2 -A 2 )+^ r (k a ,k a + k).(20) 

Using this notation we can take the w-dependence to the right hand side of Eq. ()16|) . 



i.e. 



= (v2 I 7rr 2 {^ (2) (k Q - k 6 ) (21) 

a;-^ (k 2 ,A 2 ,/i(k 2 ,k^)) l 

+ J d 2 k t r (ka, k a + k, A 2 , fi (k 2 , k 2 )) f u (k a + k, k 6 )| , 
and iterate this expression. This procedure leads to 

f (k k) - 6i2) (k " " kb) (22) 
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+ J d 2 k, 



x 



+ J rf 2 kx 



£ r .(k a; k a + k 1 ,A 2 ,/i(kg,k fc 2 )) 

(k. + kx-k,) 

uj - u ((k a + ki) 2 , A 2 , [i ((k a + ki) 2 , k; 

IC r (k a ,k a + k 1 ,A 2 ,/i(k 2 ,k 2 )) 



>t fc r (ka + k X , k a + kx + k 2 , A 2 , ((ka + k : ) 2 , k 2 



u-ljo (k 2 ,A 2 ,/i(k 2 ,k 2 )) 
A 

x / cf k 2 
x 



w - u ({ka + ki) ,A 2 ,/i((k a + k!) ,k 2 
5( 2 ) (k a + kx + k 2 - k 6 ) 



a; — ujq 



((ka + k x + k 2 ) 2 , A 2 , /i ((k a + ki + k 2 ) 2 , k 



+ 



As we want to obtain the final solution in energy space we perform the inverse 
Mellin transform of Eq. (|22|) as defined in Eq. (J2J). In this way it is possible to obtain 
the following expression for the solution in the NLLA 1 : 



/(k a ,k 6 ,A) = exp(cuo(k 2 ,A 2 ,/i(k 2 ,k 2 ))A){^(k a -k 



(23) 



+ EIIM 2 k 
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6(k 2 -\ 2 ] 
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1=0 
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i-l 



i-l 



-w ka + ^kJ ,A 2 ,/i k a + ^k, 



1=0 



1=1 



1=0 



8® K> + k a -k 6 



\Z=1 



The expression (J23"|) is only weakly A-dependent in the A — > limit. We can show 
this point by applying it to a given smooth function, $(k a ,k 6 ), and selecting those 
terms proportional to A, i.e. 

jd 2 k fe /(k a ,k fe ,A)$(k a ,k ) = $(k a ,k a ) + A|w (k 2 ,A 2 )$(k a ,k a ) (24) 
. /£(k 2 -A 2 ^ ^ ^ 



+ d'k x 



7lk 2 



aki)+ ic r (k a ,k a + ko $ (k a ,k a + ki) } + 



The A-dependence of this expression can be studied in a similar manner to the one 
applied for Eq.(18). 

1 Using the notation yo = A. 
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As a final remark we would like to point out that our solution in the LLA limit 
is consistent with similar results in the LLA in the literature jH]. If we take the 
leading-logarithmic limit of our expressions we have 

a 2 

u; (q 2 ,A 2 ) = -a s ln^, £ = a s , 77 = 0, K r (q, q') = 0. (25) 
In this case the solution takes the simple form 

/(k ,k b! A) = f^V i^ (2) (k a -k b ) (26) 



& 2 



~ JL r 2l #(k 2 -A 2 ) p-i 



n=l «=1 



V(k Q + ELik/ 



x 



<J< 2 > ($>+k a -k 6 Jj. 

This expression is equivalent to that obtained in |Hj in the LLA where the regular- 
isation of the infrared divergences was performed in a different manner. The main 
result of this Letter is the extension of this iterative method to the NLLA. This has 
been possible by using the NLLA BFKL equation in 4 + 2e dimensions to introduce 
a cut-off in the phase space. This renders the e — > limit finite and thus allows us 
to iterate the BFKL equation to obtain the solution in Eq. 



4 Conclusions 



We have presented a procedure to solve the BFKL equation for forward scattering in 
the next-to-leading logarithmic approximation. We have shown how it is important 
to use the kernel in dimensional regularisation and introduce a cut-off, A, in the phase 
space, Eq. (JTHJ). This allows us to write the solution in a compact form, Eq. (J2BJ), 
suitable for numerical studies, which will be presented in a future work. We have also 
shown the mechanism by which the solution is weakly dependent on the cut-off for 
small values of A. We would like to point out that we keep the full angular information 
in our solution by solving the equation for any conformal spin, i.e. we do not use the 
angular-averaged kernel (see Eq. (fT7])). This will allow the study of spin-dependent 
observables in the NLLA. 

Work is in progress to understand the BFKL resummed gluon Green's function 
using this approach, and to quantify the effect of those terms related to the running 
of the coupling j^j compared to the scale invariant ones. A study of the solution in 



7 



the NLLA of the N=4 supersymmetric case [TU], where the coupling does not run, is 
also under consideration. 

The ultimate goal in the application of our solution is the calculation of cross- 
sections in the NLLA for those processes where the BFKL resummation should be 
relevant. The presented solution has the advantage of not operating in the Mellin 
space, i.e. we can use the NLLA impact factors ^T] directly in transverse momentum 
space in order to make phenomenological predictions. It will, in principle, be possible 
to disentangle the structure of the final state allowing the study of, e.g., multiplicities, 
extending the work of ^2] to the NLLA. As a final point we would like to indicate 
that this approach is also suitable for studying the solution of the BFKL equation in 
the NLLA for the non-forward case when the kernel for this case has been completed. 
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